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Abstract: We investigate numerically chiral symmetry restoration at finite temperature 
in the planar limit in the deconfined phase, both when it is stable and when the system is 
supercooled. We find chiral symmetry restoration at Ty. = T^, where is the temperature 
of the deconfinement transition in pure gauge theory and T^. < in the supercooled 
deconfined phase. In the stable case the spectrum of the Dirac operator opens a gap in a 
discontinuous manner and in the supercooled phase the gap seems to vanish continuously. 
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1. Introduction. 

At finite temperature QCD undergoes qualitative changes of great physical interest. Al- 
though much is known, the complicated strongly coupled aspects of the transition region 
are not under full theoretical control. 

At a large number of colors, in the 't Hooft limit, the system simplifies somewhat Q. 
The major effects occur now in the pure gauge sector, while the fermions only react to 
these effects, without influencing them by back reaction (except, as explained later, by 
aligning the pure gauge vacuum). The purpose of this paper is to present numerical results 
on chiral symmetry restoration at finite temperature in the planar limit. Our study is at 
zero quark mass. Preliminary results have been presented last year |||]. 

At infinite Nc the free energy is temperature (T) independent at order for T < T^, 
where is the deconfinement temperature. Chiral symmetry is spontaneously broken 
and the condensate -^{'4)'4}) is nonzero and temperature independent |^, The chiral 
symmetry breakdown is reflected by a condensation of eigenvalues of the Euclidean Dirac 
operator near zero. This condensation emerges naturally in a random matrix context. 
Because of the infinite number of colors and the lack of relevance of the size of the system 
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due to large Nc reduction, one can think of the Euclidean Dirac operator (D) as a large 
random anti-hermitian matrix, whose structure is restricted only by chiral symmetry. 

, C\ 

D= 1 1.1 



-Ct 

In the spirit of Wigner's approach to complex nuclei, one is lead to write down the simplest 
probability distribution for the matrix C whose linear dimension is proportional to N^. 

P{C) oc e-«<ii-(^)^-^^^ (1.2) 

Chiral symmetry breaking is then an immediate result, giving it the appearance of a generic 
phenomenon. 

As the temperature is raised, at T = T^, a first order deconfinement transition occurs 
at all Nc > 3; Td has a finite large Nc limit ^. For high temperatures, T » T^, one 
expects chiral symmetry to be restored and, consequently, the random matrix viewpoint 
that worked for T < to become invalid. The simplest way in which chiral symmetry 
can get restored is for the Euclidean Dirac operator to open a gap at zero. One might have 
thought that this effect can be incorporated into an extended random matrix model 
but, a numerical investigation to be described later on, indicates that the types of random 
matrix models one would naturally guess will not work when chiral symmetry is unbroken. 
For T < Td random matrix theory applies also at finite iVc, but, without going to the 
planar limit the argument does not extend to high temperatures P], where there is no 



energy regime dominated by Goldstone particles |1C]. We see that going to infinite Nc does 
not help in this respect. 



2. Large Nc in the deconfined phase. 

At T > Td, in the deconfined phase, the free energy of the gluons starts depending on T. 
Feynman diagrams containing fermion loops are still suppressed by one power of so 
long as the number of flavors is kept fixed, as we do. 

In the Euclidean path integral formulation, physical finite temperature is reflected 
by the "time" direction being compactified to a circle of radius ^ and bosons/fermions 
having periodic/antiperiodic boundary conditions in the time direction. For T < Tc the 
boundary conditions are irrelevant, since the preservation of the related Z[Nc) U{1) 
global symmetry washes them out. When T > T^, the trace of parallel transport round 
the temporal circle (Polyakov loop) acquires a fixed phase and breaks spontaneously the 
associated Z{Nc). Which phase is picked is arbitrary in the absence of fermions, as all 
phases have the same gluonic energy. When fermions are present, although in general their 
contribution to the free energy is subleading, they fix the phase of the Polyakov loop to 
one specific value because the fermions break the Z(Nc) explicitly and align the vacuum. 
It is physically plausible, and supported by our numerical work, that the phase is chosen 
to make the Polyakov loop positive. This preserves CP, but also making the Polyakov loop 
negative would have. Other than aligning the vacuum, fermions have no impact on the 
distribution of the gluonic fields at infinite Nc- 
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3. Lattice setup. 

We work on a hypercubic lattice of shape L^L^. The gauge action is of single plaquette 
type. 

^=7^ E Tr[U,A^) + Ul,{x)] (3.1) 



c 



Uf^A^) = U^{x)Uu{x + fi)Ul{x + i')Ul{x) (3.2) 

We define h = = -^Tjq- = \ and take the large Nc limit with b held fixed. As usual, 
h determines the lattice spacing a and A is the 't Hooft coupling. The gauge fields are 
periodic, x is a four component integer vector labeling the site, and /i labels a direction; a 
unit vector in the fj, direction is denoted by jl. The link matrices U^{x) are in SU{N(.)- 
There is a Z^{Nc) symmetry under which 

?7^(x) ^e^C/^(x) (3.3) 

for all x with = c^. The integers are fixed, and the integers label the elements 
of the n-th. Z{Nc) group; = 0, 1, .., — 1. Changing the c^'s amounts to a local gauge 
transformation. We have Li = L2 = L3 = L. 

The Polyakov loop matrix is denoted by P4,{x) and defined by: 



P^ix) = Ui{x)Ui{x + l)U^{x + 2 • 4)..C/^(x + (L4 - 1) 



•4) (3.4) 



Under the Z{Nc) factor associated with the time direction, P^i^x) gets multiplied by a 
phase. The gauge invariant content of ^4(0;) is its set of eigenvalues (the spectrum) e * , i = 
1, 2..., Nc- The ordering is not gauge invariant, and there is a constraint that det P4{x) = 1. 
Under the Z{Nc) transformation, the set of eigenvalues is circularly shifted by a fixed 
amount. The spectrum of P4{x) and of P4{x+j - A) are the same for all j = 0, 1, 2...., L4 — 1. 

For a fixed 6 in a certain range, if L is large enough, the global Z'^{Nc) symmetry 
associated with the spatial directions is unbroken. In practice, it is even possible to work 
in metastable phases, as long as these Z{Nf^^s are maintained. 

Depending on L4 < L the Z{Nc) in the time direction may be broken or not. Alter- 
natively, for Li < L, one can break the iime,- Z{Nc) by increasing b (for L4 = L, we view 
as time the one particular, but randomly selected, direction that breaks its Z{Nc)). The 
breaking point is the deconfinement transition. There is no dependence on L on either side 
of this transition so long as b is in the range which preserves the spatial Z'^{Nc) 111]; this 
will always be the case. 

In lattice units we have aT^ = jj^: where jll] 



be{b) e{b) = ^{TrU^A^)) (3.5) 



) 51 
121 2ATT^b, 
e-^ (3.6) 
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In 1 11 1 it was found that planar gauge theory on the torus can exist in five phases, 0c,lc,..,4c; 
the deconfined phase is the Ic phase in this notation. 

The fermion action is given by the overlap Dirac operator which preserves chiral sym- 
metry exactly. This choice makes it possible to pose the question of spontaneous chiral 
symmetry breaking in a clean way. 

The massless overlap Dirac operator [12|, Do, is defined by: 



— 2 

F-1 = = ^^y^^ = sign(i/^(M))75 (3.7) 

Hw{M) is the Wilson Dirac operator at mass M, which we shall choose as M = —1.5. M 
should not be confused with a bare quark mass. 



Hy,{M) = 75 



(3.8) 



The matrices are the lattice generators of parallel transport and depend parametrically 
and analytically on the lattice links U^{x). 

The internal fermion-line propagator, is not needed at infinite Nc, since fermion 
loops are suppressed. For fermion lines continuing external fermion sources we are allowed 
to use a slightly different quark propagator [O] defined by: 



I l-V , , 

A = TT7 

A = —A'^ and anticommutes with 75. The spectrum of A is unbounded, but is determined 
by the spectrum of V which is restricted to the unit circle. Up to a dimensionful unit, A 
should be thought of as a lattice realization of the continuum massless Dirac operator, D: 

2\M\A ^D = 7^a^ + (3.10) 

Our main observable will be the smallest eigenvalue of the non-negative matrix —A^, 
which is the discrete version of D^D, where D is the continuum, Euclidean, Dirac operator 
in a fixed gauge background. The gauge background varies according to the pure gauge 
action and we shall look at the induced probability distribution of the smallest eigenvalues 
of -A\ 

Let Ai^2 be the two lowest eigenvalues of V—A'^: The dimensionless gap as a function 
of the dimensionless temperature t, g(t), is defined as the average over gauge configurations 
of Xi{h, L,Li)Lc{h). The dimensionless temperature itself is defined as Lc{b)/Li. 

At infinite Nc, g will vanish when the symmetry is spontaneously broken. If we find 
that g is nonzero, we know that chiral symmetry is restored. This is true because the 
single way chiral symmetry transformations can avoid inducing naive relations among cor- 
relation functions of physical observables is by the matrix A~^ becoming singular with finite 
probability. If g > almost surely, this cannot happen and chiral symmetry is preserved. 

As we vary L4 or b, g will change. If the change is discontinuous to or from zero, 
the transition is of first order; otherwise, it is continuous. We shall find that g jumps 
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discontinuously from zero to a finite value exactly when the pure gauge field undergoes the 
phase transition. Also, we shall see some evidence that when we lower the temperature 
T below T^, but stay in a supercooled deconfined phase, eventually, g goes to zero at 
some temperature < T^^; this chiral symmetry breaking transition, occurring as the 
temperature is lowered in the supercooled deconfined phase, seems continuous within our 
numerical resolution. 

4. Vacuum alignment. 

We wish first to check that indeed the vacuum of the gauge fields is selected by aligning 
with the fermions in the manner discussed earlier. 

Numerically, we would like to show that the (positive) fermion determinant is indeed 
favoring a positive Polyakov loop in the deconfined phase. We do this as follows: 

We let the lattice system pick an arbitrary phase for its Polyakov loop in the deconfined 
phase. We now define the fermions with twisted boundary conditions relative to this 
phase in the time direction. That is, the fermions obey the boundary condition "0(0) = 
—il)(^-^)e'^^^~'^\ where the Polyakov loop has phase e*"^. We now intend to show that the 
fermion determinant is maximal when = 0. 

A complete computation of the determinant is too expensive and might be an overkill. 
We accept the hypothesis that the determinant is maximal when the gap g is. After all, 
the eigenvalues of the Dirac operator repel and a larger g simply would push all eigenvalues 
to slightly higher values, hence increasing the determinant itself. 

In summary we end up focusing on the gap 5 as a function of the angle 6. As expected, 
we obtain a periodic function of 9 with a local maximum at = and local minima at 
6 = zbTT, symmetric under 6 —9 and monotonically decreasing from its value at = to 
its value at = vr. 

CP invariance implies ^ = at = 0, vr; figure || also shows that the region of 
maximum is connected by an approximately linear segment to the regions of the minima. 

This linearity can be understood if one accepts a static approximation, which is plausi- 
ble for high enough temperatures. In this static approximation the gauge field in the time 
direction is taken as a constant, and the gauge fields in the space directions are taken as 
time independent. In the continuum, the spectrum of the operator 74D (where 4 labels the 
time direction) has its spectrum linearly shifted by assuming now that this shift gets 
transmitted almost intact to the lower bound of the operator D, we obtain the linearity 
of the gap g with a predicted slope. 

Note however, that as the temperature is lowered towards T^, T gets replaced by a 
smaller temperature, T^ff < T. Still, the gap is maximal at = and vanishes at = ivr, 
and the linear portion is always present. 

One might have expected to find the gap vanishing before 9 reaches itvr, for the follow- 
ing reason: Parallel transport round the time direction results in a unitary matrix whose 
spectrum has support on a connected arc on the unit circle, centered at unity. Because 
of the finite width, one could have imagined that in some color direction fermions are ef- 
fectively defined with a twisted boundary condition and the static approximation would 
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Figure 1: The gap as a function of twist. 



have predicted an earlier point where the gap associated with that color direction vanishes. 
We do not see this happening, indicating that there is no sense in which the eigenvalues 
associated with the Polyakov loop can be used to coherently select special orientations in 
color space for the fermions. All that one can expect is an incoherent effect, in which the 
spread in eigenvalues influences averages over all colors - a point we shall return to later. 

In the static approximation was used to motivate a variation on the random matrix 
model describing the Dirac operator in the confined phase. The Dirac operator is now split 
into blocks labeled by the Matsubara integer n, and each block has a structure 



/ C + z[i2n + l)n-e]T\ 

" " V-Ct + z[(2n + lK-0]T J ^ ' 

with a common matrix C distributed according to equation pT^ . 

In this picture the Polyakov loop is taken as a unit matrix with one overall phase. If 
we generalized the above random model by further splitting into blocks associated with 
each color we would obtain a prediction that the g{9) ought to vanish before 9 reaches ivr, 
but this is somewhat implausible and indeed does no happen as pointed out above. 
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5. A possible implication of the spread in the eigenvalues associated with 
the Polyakov loop. 

We present here a suggestion that the spread in eigenvalues explains why one obtains an 
effective temperature T^ff < T when one gets closer to from above. 

Consider a free Dirac fermion in continuum with twisted boundary condition 9. Let 
/ be the contribution of this one fermion to the free energy density in Euclidean space. 
Standard manipulations [0] produce 

(5.1) 

For 9 = we recover the well known expression. As 9 increases toward vr the coefficient 
of decreases and eventually even changes sign. Indeed, at = vr we have boundary 
conditions appropriate for a boson, but since we kept the determinant at a positive power, 
we used for it Fermi statistics. Had we used the right statistics, the determinant would 
have been at a negative power and then the contribution to / would have had the normal 
sign. 

We now speculate that the fact that parallel transport round the compact direction 
is best described with a phase drawn from a distribution symmetric about zero but not 
of zero width, effectively implies an averaging over 9 in the above equation. In turn, this 
could be viewed as an effective reduction of the temperature, if one insists on keeping the 
classical value for the prefactor. 

For large T, T >> T^, the width of the distribution shrinks to zero and the effect 
goes away. However, for T close enough to T^, this indicates that / ~ c{T)T^ where c{T) 
decreases from its classical value at T = cxd slowly. The fermionic contribution to the free 
energy would therefore appear almost noninteracting, however with a coefficient that is 
slightly off. 

Much more is needed to see if this speculation bears out. 

6. Supercooling. 

So long as we have confinement we ought to have spontaneous symmetry breakdown at 
infinite Nc for theoretical reasons We also know from numerical work that the decon- 
finement transition in planar QCD is of first order. 

We now ask what happens to chiral symmetry as we supercool the deconfined phase, to 
temperatures T < T^. The gauge coupling should increase, and eventually, chiral symmetry 
could break again, this time without confinement. We address this question numerically. 
Because of the need to extrapolate into the metastable phase, the conclusion is quite 
tentative. 

We see an indication that there is a chiral symmetry breaking temperature Ty. < T^, 
where chiral symmetry is broken in the supercooled deconfined phase. Our numerical 
finding is consistent with a continuous transition at T^. 
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Our numerical results are collected in table ||. All the results are in the Ic phase and we 
have used anti-periodic boundary conditions for fermions with respect to the Polyakov loop 
in the broken direction. We studied five different couplings, namely, b = 0.35, b = 0.3525, 
b = 0.355, b = 0.3575 and b = 0.36. We will use the following central value for the critical 
sizes at these couplings: 

Lc(0.35) = 5.97; ^^(0.3525) = 6.46; ^^(0.355) = 6.96; 
Lc(0.3575) = 7.47; Lc(0.36) = 8.01. (6.1) 

We performed a careful study at 6 = 0.3525 and L4 = L = 6 and convinced ourselves that 
the large N limit is obtained for values listed in the table. With the exception of the entry 
at 6 = 0.355, N = 23, L = 8 and -L4 = 7, all entries are definitely in the stable Ic phase. 
It is possible that the b = 0.355, = 23, L = 8 and L4 = 7 entry is in the supercooled Ic 
phase. 

The data at the largest A'' from table || was used to obtain the dimensionless gap, g{t), 
as a function of the dimensionless temperature t. Figure || shows all the data for t < 2.8. 
There is some spread of points, indicative of order lattice corrections, but the data seems 
to condense toward a line to which we assign a physical meaning in the continuum. We 
show a two parameter fit to 

g = 1.76Vt - 0.93 (6.2) 

The square root behavior was imposed on the fit. 

This indicates a continuous chiral restoring transition in the supercooled phase at T ~ 
0.93rrf. We guess the transition ought to be continuous because the density of eigenvalues 
of the Dirac operator at zero determines it, and so long as we stay in the deconfined phase 
(metastable regions included) the dependence on the temperature of the gauge background 
has no a priori reason to change discontinuously; one could view the temperature as entering 
only through the gauge background, but this point might be debated. 

7. No naive random matrix description. 

Intuitively, one could argue that large Nc alone, without the additional input from the 
viability of an effective chiral Lagrangian description is enough to motivate a random 
matrix description of spectral properties of the Dirac operator in planar QCD. If this were 
true, one ought to be able to describe the spectrum of the Dirac operator by some random 
matrix model even at temperatures where chiral symmetry is restored. 

We tried to see if this is possible by looking at the correlation between the fluctuations 
of the lowest eigenvalues of the Dirac operator. These fluctuations are about the mean, 
so independent of the gap g, which is the main difference between the case where chiral 
symmetry is broken and we know that random matrix theory works and where we are when 
chiral symmetry is restored. 

If indeed there is some random matrix model of the usual type considered the cor- 
relation between such fluctuations is mainly governed by level repulsion. We have two 
examples of such random matrix models: one where chiral symmetry is broken and the 
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Table 1: Results for the lowest two eigenvalues and their correlation c as defined in ( |7.l[ ) are 
listed for several values of 6, i, L/^ and N in the Ic phase. 



other where it is not, due to a large enough shift of C (which occurs for T large enough) 



in the formula for D in equation 4.1 



We estimate by Monte Carlo methods the correlation, c between the two lowest eigen- 
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Figure 2: The gap in units of Td as a function of temperature in units of T^. 



values: 



((Ai-(Ai))(A2-(A2))) 



(7.1) 



(Al))2)((A2-(A2))2) 

In the two random matrix models mentioned before c is between ^ (deep in the symmetric 
phase) and ^ (broken phase), c is bounded from above by 1 by a simple positivity argument. 
The numbers in table [l| show that c is quite close to unity deep in the symmetric phase and 
and, although it drops down a little as the temperature decreases, c remains well above 0.5 
even close to i = 1. This indicates a much stronger correlation between the fluctuations of 
the lowest eigenvalues of the Dirac operator than one would expect in any random matrix 
model. The correlation we find is so strong as to imply an almost rigid relationship between 
the two lowest eigenvalues. 

If there exists a random matrix model that applies to the case when chiral symmetry is 
restored, we are missing an essential ingredient, which we speculate might have something 
to do with the spectra of the eigenvalues associated with the Polyakov loops. 



8. Comparison to holographic models. 

Recently several papers discussing finite temperature features of models that bear various 
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degrees of semblance to QCD have appeared ||16|, |17|, These models are in the contin- 
uum, but admit dual descriptions which allow control of the strong coupling regime in the 
planar limit. 

There are several models involving a set of branes which produce the gauge fields 
and the dual gravitational background and "probe" branes, that have no effect on the 
background, which produce the quark fields. One finds deconfining and chiral symmetry 
breaking transitions, all of first order. 

In particular, |18| addresses the issue of a chiral symmetry restoration transition in the 
supercooled deconfined phase and finds it is of first order. It seems that a crucial ingredient 
is the presence of an extra scale, beyond Aqcd, which allows for a more complex phase 
diagram. Nevertheless, it is noteworthy that some of these solvable situations are members 
of continuously varying families, which also contain ordinary QCD, albeit in a regime that 
is not under control in the dual variables. 

Apart from the order of the transition in the supercooled phase, which admittedly is 
only tenuously determined here, there is general agreement between our findings here and 
the holographic results. 



9. Conclusions 



Our main conclusion is unsurprising: at infinite Nc the strong first order deconfinement 
transition induces also immediate chiral symmetry restoration. 

At a more minor level we obtained some less expected results: If we supercooled the 
deconfined phase chiral symmetry would still break spontaneously and seemingly does so in 
a continuous transition. Moreover, that transition would occur at a temperature which is 
only seven percent below the deconfinement transition. It is well known that such "pseudo- 
transitions" are indicative of complex dynamics even in the stable phase. We also saw that 
there is something fundamentally different in the statistics of the eigenvalues between the 
Dirac operator and typical random matrix theory models. It would be worthwhile to 
understand why this is so and what it means. In addition we were led to a heuristic 
explanation for why the coefficient c in the free gas law for the free energy, / = cT^, is 
weakly temperature dependent for temperatures close to and smaller than the classical 
value. The explanation viewed this as a consequence of the fact that the eigenvalues 
associated with the Polyakov loop are distributed over an arc of a certain width, and only 
at infinite temperature does the Polyakov loop become the trace of a matrix proportional 
to the unit matrix. 

The new tools of holographic duals provide a wealth of exactly soluble strongly coupled 
theories where similar phenomena occur and it would be useful to find observables and 
models where both lattice and holographic methods apply simultaneously. 
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